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^ ' In this work we consider the influence of mass media in the dynamics of the 

^ . two-dimensional Sznajd model. This influence acts as an external field, and it is 

O ' introduced in the model by means of a probability p of the agents to follow the media 

^ ■ opinion. We performed Monte Carlo simulations on square lattices with different 

sizes, and our numerical results suggest a change on the critical behavior of the 
Q^l model, with the absence of the usual phase transition for p >~ 0.18. Another effect 

of the probability p is to decrease the average relaxation times r, that are log- 
Ps| . normally distributed, as in the standard model. In addition, the r values depend 

I on the lattice size L in a power-law form, r ~ L", where the power-law exponent 

ly-^ ' depends on the probability p. 
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' 1 Introduction 



Social dynamics have been studied through statistical physics techniques 
in the last twenty years. Among the studied problems, we can cite models of 
cultural P, language [2] and opinion dynamics [3P|5] (for a recent review, see 
[6]). These kinds of models are interesting to physicists because they present 
order-disorder transitions, scaling and universality, among other typical fea- 
tures of physical systems [6]. 
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One of the most studied models of opinion dynamics in the last years is the 
Sznajd model [5|7] . The original Sznajd model [5] consists of a chain of sites 
with periodic boundary conditions where each site (individual opinion) could 
have two possible states (opinions) represented in the model by Ising spins 
("yes" or "no"). A pair of parallel spins on sites i and i + 1 forces its two 
neighbors, i — 1 and i + 2, to have the same orientation (opinion), while for an 
antiparallel pair [i, i + the left neighbor {i — 1) takes the opinion of the spin 
i + 1 and the right neighbor [i + 2) takes the opinion of the spin i. In this first 
formulation of the model two types of steady states are always reached: com- 
plete consensus (ferromagnetic state) or stalemate (anti- ferromagnetic state), 
in which every site has an opinion that is different from the opinion of its 
neighbors. However, the transient displays an interesting behavior, as pointed 
by Stauffer et al. [8j. 



A more interesting situation was studied in [S], where the model was defined 
on a L X L square lattice. The authors in |8] considered not a pair of neighbors, 
but a 2 X 2 plaquette with four neighbors. Considering that each plaquette 
convince its eight neighbors if all four center spins are parallel, and that the 
initial density of up spins is d = 0.5, the authors found that the system reaches 
the fixed points with all up or all down spins with equal probability. For d < 0.5 
(> 0.5) the system goes to a ferromagnetic state with all spins down (up) in 
all samples, which characterizes a phase transition at d = 0.5 in the limit 
of large L. This phase transition separates two distinct states of the system: 
for d < 0.5 the system never reaches full-consensus states with all spins up, 
whereas for d > 0.5 the consensus is always reached. Other formulations of the 
two-dimensional model, considering for example memory [9], reputation |10j . 
diffusion of agents [llj or a random dilution of the lattice |12], also present 
similar phase transitions. 

The utility of the Sznajd model goes beyond the basic description of a commu- 
nity [7]. It was applied, for example, to politics. In 1999, Costa-Filho et. al. [13] 
showed that distribution of votes per candidate for the 1998 elections in Brazil 
follow a power law with exponent ~ —1.0. Based on the Sznajd model, two 
models with more than two opinions were proposed, considering the dynamics 
on square and cubic regular lattices, and also in complex Barabasi networks 
[2]. The authors successfully reproduced the distribution of the number of 
candidates according the number of votes they received in Brazilian elections 
[T5|16|17j . The Sznajd model was also applied to marketing, where it was 
considered advertising and a competition of two different products, and to 
finance, where the authors found a good agreement between some characteris- 
tics of the price trajectories (like returns, for example) and the simulations [7j. 
It was also claimed that some other characteristic of Sznajd-like models may 
also be present in real social systems, like the power-law relationship between 
the time need to reach the fixed point (the complete consensus), called the 
relaxation time, and the system size |llj . 
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However, we can observe that the above-mentioned works did not take into 
account the effects of an external field in this phase transition. In a real com- 
munity, external effects may be, for example, the mass-media influence. It is 
well-known that, in real life, the mass media (television, radio, ...) has a great 
influence in the population, and people tend to keep or change their opinions 
on any question according to that influence. These effects were considered 
in some social models, specially in the Axelrod model of cultural diversity 
|18|19|20|21|22|23] . and interesting results were found. Thus, in this work we 
consider the effects of mass media in the opinion dynamics of the Sznajd so- 
ciophysics model. This influence is introduced in the model by means of a 
probability p of the agents to follow the media opinion. We performed Monte 
Carlo (MC) simulations for different lattice sizes, and our results suggest that 
the system undergoes a phase transition, as well as in the original Sznajd 
model, for values of the probability p <~ 0.18, and that the relaxation times 
are log-normally distributed. In addition, the average relaxation times r de- 
pend on the lattice size L in a power- law form r ~ L", where a is a function 
of p. 

This work is organized as follows. In Section 2 we present the microscopic 
rules that define the model. The numerical results are discussed in Section 3, 
and our conclusions are presented in Section 4. 



2 Model 

We have considered the Sznajd model defined on a square lattice with linear 
size L and periodic boundary conditions. The lattice sites were numbered by 
one index i, i = 1,2,...,N, where = is the total number of agents in 
the population. We assign an Ising variable to each site, Sj = ±1, representing 
the two possible opinions of each agent. At each time step, the following 3 
microscopic rules control our model (see the schematic representation of Fig. 
CD: 

(1) We randomly choose a 2 x 2 plaquette of four neighbors. 

(2) If all four center spins are parallel, the eight nearest neighbors are per- 
suaded to follow the plaquette orientation. 

(3) If not all four center spins are parallel, we consider the influence of mass 
media: each one of the eight neighbors follows, independently of the oth- 
ers, the media opinion with probability p. 

Notice that we update the agents' states in a random sequential order (asyn- 
chronous updating). After the above three steps, we count one MC step in 
our model. In addition, we will consider that the mass media is favorable to 
the opinion i.e., in the above-mentioned rule (3), each neighbor (inde- 
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Fig. 1. Schematic representation of a plaquette (spins inside the rectangle) and its 
eight neighbors (full-hne spins). Notice that the other agents (dotted- hne spins) do 
not participate in the dynamics. 

pendently of the others) will change his opinion to +1 with probability p if 
he was not persuaded by the plaquette agents (i.e., if not all four plaquettes' 
spins are parallel). In other words, we will take into account that the influence 
of a group of agents is stronger than the influence of media. In fact, we can 
imagine that in the real world people tend to be more influenced by friends, 
relatives and colleagues (represented in the model by the plaquettes' agents), 
among others, in comparison with the influence of the media (represented in 
the model by the probability p). Notice that for p = we recover the standard 
Sznajd model defined on the square lattice [S]. 

As pointed by Stauffer in [21], we can imagine that each agent in the Sznajd 
model carries an opinion, that can either be +1 or —1 and that represents 
one of two possible opinions on any question. For example, in an election, the 
opinion +1 may represent the intention of an individual to choose a certain 
candidate A, whereas the opinion —1 may represent the intention to choose 
another candidate B. Thus, the inclusion of the probability p in the model will 
take into account the influence of TV and/or advertisement, for example, to 
provoke a change in the choice (i.e., the candidate) of that individual. 



3 Monte Carlo Simulation 

Following the previous works on the Sznajd model, we can start studying the 
time evolution of the magnetization per site. 



where = is the total number of agents in the population and = ±1. 
In the standard Sznajd model defined on the square lattice [8], if we consider 
an initial density of up spins d = 0.5 and that a 2 x 2 plaquette with all spins 



m = — 




(1) 
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Fig. 2. (Color online) Time evolution of the magnetization per site m, Eq. ([1]), for 
L = 73 and different samples. The parameters are p = 0.03, d = 0.38 (a) and 
p = 0.1, d = 0.13 (b). Observe that in both figures the consensus with all spins 
pointing up is achieved in some samples, even for d « 0.5, that is the critical 
density for the standard model [8j. 

parallel persuades its eight neighbors, the system reaches the fixed points with 
all up or all down spins with equal probability. For d < 0.5 (> 0.5) the system 
goes to a ferromagnetic state with all spins down (up) in all samples, which 
characterizes a phase transition at a critical value dc = 0.5 in the limit of 
large L. We can observe that for p > this result changes (see Fig. [2]). In 
Fig. M^a.) we exhibit the simulation data for p = 0.03 and initial density of 
up spins d = 0.38. In this case, we observe that in some realizations of the 
dynamics the system reaches consensus with all spins up, even for d < 0.5. In 
the case of a greater value of p, for example p = 0.1, the difference between our 
model and the standard one [8] is more pronounced, and a strong effect of the 
media influence is observed: even for a very small value of the initial density 
of up spins {d = 0.13) we observe that the full-consensus state with m = 1 
is obtained in some samples [see Fig. El^b)]. Notice that as in the original 2D 
Sznajd model [8], the system always reaches consensus states with all up or 
all down spins, i.e., there are no states with magnetization per spin \m\ < 1. 
However, the presence of an external field (mass media) changes the critical 
behavior of the 2D Sznajd model. In addition, if a phase transition also occurs 
in our case, the position of the transition points dc will depend on p. 

The next step is to analyze if the phase transition observed in the standard 
2D Sznajd model also occurs in our version of the model. In other words, 
given a certain value of the probability p, we want to determine if there is a 
critical value d^ above which the system always reaches full-consensus states 
with m = 1. For this purpose, we have simulated the system for different 
lattice sizes L and typical values of the probability p, and we have measured 
the fraction of samples which show all spins up when the initial density of up 
spins d is varied in the range 0.0 < d < 1.0. In other words, this quantity 
/ gives us the probabihty that the population reaches consensus for a given 
value of d. The results for the largest size considered, L = 121, are exhibited in 
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Fig. 3. (Color online) Fraction / of samples which show all spins up when the initial 
density of up spins d is varied in the range 0.0 < d < 1.0 for L = 121 and typical 
values of the probability p. In the inset we show results for values in the range 
0.12 < p < 0.2. Observe that the transition does not occur for p = 0.2, i.e., we have 
/(d = 0) = and /(d > 0) = 1. We have considered 300 samples for each point in 
the graphics. 

Fig. [3l where we have considered 300 samples for each point. We can observe 
the transition from / = to / = 1 at different values of d, for a given value of 
p. However, for p = 0.2 the transition does not occur, i.e., we have that / = 1 
even for very small values of the initial density of up spins like d = 0.01 (see 
the inset of Fig. [3]). The exact threshold value pthres above which the system 
does not present the phase transition is difficult to determine numerically. 
Nonetheless, our simulations suggest that the phase transition only occurs in 
our formulation of the Sznajd model for p <~ 0.18. 

To see this picture in more details, we performed simulations for L = 121 
and d = 0.01, and the results obtained were / = 0.97 for p = 0.18, / ~ 0.98 
for p = 0.19 and / = 1.0 for p = 0.2, whereas for p = 0.16 and p = 0.15 
we obtained / = 0.69 and / = 0.46, respectively. We can obtain the same 
conclusion monitoring the time evolution of the magnetization per site m. 
For example, in Fig. H] we exhibit the magnetization m as a function of the 
MC time for d = 0.01, p = 0.2 and different samples. Observe that even 
for a very small value of the initial density of up spins the full consensus is 
achieved in all samples, suggesting the absence of a phase transition (similar 
results were found for 0.18 and 0.19). In other words, the mass media (i.e., 
the probability p) acts in the system in a way that it changes the opinion of 
the majority of the agents in the population, and the full-consensus steady 
states with m = 1 are always achieved for p >~ 0.18. As above-discussed, 
we are interested in the transition between a phase where the full-consensus 
states with all spins pointing up never occur and a phase where this consensus 
always occurs. In this sense, the main effect of the mass media is to eliminate 
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Fig. 4. (Color online) Time evolution of the magnetization per site m for L = 73, 
p = 0.2, d = 0.01 and different samples. Observe that for this value of the probability 
p, even for d = 1% the full consensus with all spins pointing up is achieved in all 
samples, suggesting the absence of a phase transition. 



Fig. 5. (Color online) Log- log plot of the histogram of relaxation times for L = 73, 
d = 0.5 and some values of the probability p, obtained from 10^ samples. The dis- 
tribution is compatible with a log-normal one for all values of p, which corresponds 
to the observed parabola in the log-log plot. 

the phase transition for p >~ 0.18. 

An usual quantity measured in MC simulations of Sznajd models is the re- 
laxation time, i.e., the time needed to find all the agents at the end having 
the same opinion. The distribution of the number of sweeps through the lat- 
tice, averaged over 10^ samples, needed to reach the fixed point is shown in 
Fig. |5l We can see that the relaxation times distribution is compatible with 
a log-normal one for all values of p, which corresponds to a parabola in the 
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Fig. 6. (Color online) Average relaxation times r, over 10^ samples, versus lattice 
size L in the log-log scale for d = 0.5 and some values of the probability p. The 
straight lines are fittings that present different slopes, i.e., r depends on L in the 
power-law form t ^ L^, where a is a function of p. 

log-log plot of Fig. |5l In addition, the effect of the probability p is to decrease 
the relaxation times. In fact, the peaks of the distributions in Fig. |5]move to 
lower values of the relaxation times when we increase p. Notice that relaxation 
times log- normally distributed were observed in other versions of the Sznajd 
model pilUfllf25] . 

Based on the relaxation times distribution, we measured the average relaxation 
times T (also over 10'^ samples) as a function of the lattice size L. The results 
are exhibited in Fig. |6l in the log-log scale. As in previous works on the Sznajd 
model [IO|TT] . we have found a power- law relation between these quantities. 




but with different exponents a, depending on the value of the probability p. 
Observe in Fig. [6] that the average relaxation times are smaller for increasing 
values of the probability which is a consequence of the effects of p on the 
distribution of the relaxation times. In fact, we can observe that the peaks of 
the distributions in Fig. [5] move to lower values of the relaxation times when 
we increase as discussed above. 

In Fig. [7| we exhibit the exponent a as a function of p, and we can observe 
that a decreases for increasing values of p. It is also shown (see the inset of 
Fig. [7]) the same data in the log-log scale. Fitting data, we obtained another 
power-law behavior. 
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Fig. 7. (Color online) Exponent a of the average relaxation time, Eq. ([2]), as a 
function of the probability p. In the inset we show the same data plotted in the 
log- log scale. Fitting data, we obtained a 'p~^-^^ . 

It is important to say that the results presented in this work are robust against 
changes in some details of the model. For example, if we consider plaquettes 
with two agents (instead of four) that try to convince their six neighbors 
(instead of eight), we also observe the same results. 



4 Conclusions 

In this work, we have studied a modified version of the two-dimensional Sznajd 
consensus model. In particular, we have introduced the effects of mass media 
on the opinion dynamics of the model. This influence acts as an external 
field, and it is introduced in the model by means of a probability p of the 
agents to follow, independently of the others, the media opinion. The standard 
Sznajd model is recovered for p = 0. We performed Monte Carlo simulations 
of the model defined on square lattices, and our results are summarized in the 
following. 

We have analyzed in detail the time evolution of the magnetization per spin 
for different values of the probability p. The simulations suggest that the sys- 
tem undergoes a phase transition only for values of p <~ 0.18. This result was 
confirmed by measurements of the fraction of samples which show all spins up 
when the initial density of up spins d is varied, that is the usual order param- 
eter for the Sznajd model. Thus, the main effect of the probability p, that acts 
in the model as a mass-media (external) influence, is to eliminate the phase 
transition. This behavior is similar to the effect of an external uniform mag- 
netic field H on the 2D Ising model. The difference is that for the Ising model 
any value H > provokes the destruction of the ferromagnetic-paramagnetic 
phase transition [26], whereas in our case it is necessary an "external field" 
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p ~ 0.18 to eliminate the phase transition of the 2D Sznajd model. 

Observe that the inclusion of simple ingredients modify drastically the be- 
havior of the Sznajd model. In fact, as discussed in the Introduction, the 
one-dimensional Sznajd model does not present a phase transition. However, 
the inclusion of a temperature provokes the appearance of a transition even 
in the ID case ^7\. Considering the two-dimensional model, the simplest for- 
mulation exhibits a phase transition [8j, whereas the inclusion of an external 
field destroys this transition, as discussed in this work. 

In order to measure the relaxation times of the model, we performed Monte 
Carlo simulations on square lattices with linear sizes up to L = 121 and typi- 
cally 10^ samples. As in the standard model [8], we found that the relaxation 
times are log-normally distributed. In addition, the average relaxation times 
r could be characterized by power laws, r ~ L", where a depends on p. Based 
on the numerical values of a, we found that the dependency is also a power 
law, a ~ which implies that the average relaxation times decrease for 

increasing values of p. 
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